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ABSTRACT 

We present spectral classifications for 438 B-type stars observed as part of the VLT-FLAMES Tarantula Survey (VFTS) in the 
30 Doradus region of the Large Magellanic Cloud. Radial velocities are provided for 307 apparently single stars, and for 99 targets 
with radial-velocity variations which are consistent with them being spectroscopic binaries. We investigate the spatial distribution of 
the radial velocities across the 30 Dor region, and use the results to identify candidate runaway stars. Excluding potential runaways 
and members of two older clusters in the survey region (SL639 and Hodge 301), we determine a systemic velocity for 30 Dor of 
271.6 ± 12.2 km s _1 from 273 presumed single stars. Employing a 3-cr criterion we identify nine candidate runaway stars (2.9 % of the 
single stars with radial-velocity estimates). The projected rotational velocities of the candidate runaways appear to be significantly dif¬ 
ferent to those of the full B-type sample, with a strong preference for either large (> 345 km s _1 ) or small (< 65 km s _1 ) rotational ve¬ 
locities. Of the candidate runaways. VFTS 358 (classified B0.5: V) has the largest differential radial velocity (—106.9 ± 16.2 km s _1 ), 
and a preliminary atmospheric analysis finds a significantly enriched nitrogen abundance of 12 + log (N/H) > 8.5. Combined with a 
large rotational velocity ( v e sin i = 345 ± 22 km s _1 ), this is suggestive of past binary interaction for this star. 

Key words, open clusters and associations: individual: NGC 2060. NGC 2070, Hodge 301, SL 639 - stars: early-type - stars: funda¬ 
mental parameters 


1. Introduction 

The VLT-FLAMES Tarantula Survey (VFTS) is a European 
Southern Observatory (ESO) Large Programme that has ob¬ 
tained multi-epoch optical spectroscopy of over 800 massive 
stars in the 30 Doradus region of the Large Magellanic Cloud 
(LMC). An overview of the survey, including details of the stel¬ 
lar photometry, spectroscopic observations and data reduction, 
was given by Evans et al. (2011; hereafter Paper I). 

The VFTS observations are the first optical spectroscopy for 
many of the targets, and provide a wealth of new morphologi- 
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cal information on the high-mass population of 30 Dor, which 
will underpin quantitative analyses of the spectra in future pa¬ 
pers. Detailed classifications for the 352 O-type spectra (with 
discussion of notable morphological groups) were presented by 
Walborn et al. (2014), while classifications for the smaller sam¬ 
ples of Wolf-Rayet and late-type stars were given in Paper I. 
Here we present spectral classifications for 438 B-type objects 
from the survey. Of these, only 74 (17%) have previously pub¬ 
lished classifications based on digital data, many of which were 
obtained at a lower spectral resolving power. 

To complement the spectral classifications, we also present 
estimates of radial velocities for the B-type sample. These es¬ 
timates provide insights into the different populations within 
30 Dor, and enable identification of candidate runaway stars that 
have been ejected from their birth sites. In the context of run¬ 
away stars, the size of the VFTS sample is unprecedented for a 
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single, young star-forming region. The survey has revealed ~20 
O-type runaway candidates (Evans et al. 2010, Sana et al., in 
prep.), with an apparent connection with rapid rotation (Walborn 
et al. 2014, Sana et al., in prep.). There are also a number of stars 
with moderate (projected) rotational velocities in the peripheral 
regions of 30 Dor, which could perhaps be very rapid rotators 
(but with low inclinations; see Walborn et al. 2014, for a dis¬ 
cussion of their spectra). For completeness, we note that two 
additional 02-type stars near 30 Dor, Sk—68° 137 and BI253 
(= VFTS 072), were suggested as runaways by Walborn et al. 
(2002a) on the basis of their remote locations for their apparent 
youth. 

In his pioneering study, Blaauw (1961) noted that the frac¬ 
tional incidence of Galactic runaways is larger among O stars 
than among B stars, a result substantiated by subsequent authors, 
and which also emerged in the simulations from Portegies Zwart 
(2000). In this paper we therefore complement the studies of the 
O-type runaways in the VFTS with analysis of the B-type stars. 
In particular, the number and mass spectrum of the ejected run¬ 
aways from 30 Dor will provide important constraints on models 
of cluster formation and evolution (e.g. Fujii & Portegies Zwart 
2011). The estimated radial velocities presented here, combined 
with those from the O-type stars (Sana et al. 2013), will also pro¬ 
vide an important input for comparisons with theoretical predic¬ 
tions, e.g. in the context of simulations to investigate the initial 
conditions leading to the distinct clump identified to the north¬ 
east of R136 (Sabbi et al. 2012). 

In Section 2 we summarise the relevant features of the VFTS 
data and their post-processing. Section 3 sets out the framework 
within which the spectra were classified and discusses notable 
objects. Section 4 details the methods we used to estimate radial 
velocities for each target, followed by a discussion of members 
of the two oldest clusters in the survey region (Hodge 301 and 
SF639) in Section 5. Finally, in Section 6, we use our velocity 
estimates to identify candidate single-star and binary runaways 
and discuss their nature and potential origins. In parallel to this 
study, Dunstall et al. (in prep.) have exploited the multi-epoch 
VFTS observations to investigate the multiplicity properties of 
the B-type sample. 

2. Observations and data processing 

All of the VFTS data discussed here were obtained using 
the Medusa-Giraffe mode of the Fibre Farge Array Multi- 
Element Spectrograph (FFAMES) instrument on the Very Farge 
Telescope (VET). The Medusa fibres (which subtend T/2 on the 
sky) were used to relay light from up to 130 targets simultane¬ 
ously to the Giraffe spectrograph (see Pasquini et al. 2002). 

Full details of the observational strategy and reduction of the 
data are given in Paper I. In brief, the ESO Common Pipeline 
Fibrary FFAMES recipes were used to undertake bias sub¬ 
traction, fibre location, summed extractions of each object, di¬ 
vision by a normalised flat-field, and wavelength calibration. 
Subsequent processing included correction of the spectra to the 
heliocentric frame, sky subtraction, and rejection of significant 
cosmic rays. 

The Medusa observations used three of the standard Giraffe 
settings (FR02, FR03, and HR15N), providing spectral cover¬ 
age of A3960-5071 A (at resolving powers, R, of 7 000 to 8 500), 
and A6442-6817 A (at R = 16 000). These resolutions are greater 
than for most published standards so, for the purposes of classifi¬ 
cation, the data were degraded to an effective resolving power of 
R = 4 000, ensuring consistency with the approach to the classi¬ 
fication of the O-type stars by Walborn et al. (2014). 


For stars that appear to be single (Dunstall et al., in prep.) 
the FR02 and FR03 spectra were normalised then merged, with 
the cr-clip levels for cosmic rejection and relative weightings 
based on the signal-to-noise (S/N) ratios of the individual spec¬ 
tra. Finally, for classification, the merged FR02 and FR03 spec¬ 
tra were combined, then smoothed and rebinned to R = 4 000. 
For stars identified as having significant but relatively modest ve¬ 
locity variations (A v r < 40 km s -1 ), the individual spectra were 
similarly co-added and degraded (as the variations were within 
the effective velocity resolution of the rebinned data). For single- 
lined binaries with A v r > 40 km s^ 1 , the single-epoch FR02 
and FR03 spectra with the best S/N ratio were combined and 
then smoothed/rebinned. We were careful in these cases to en¬ 
sure that any velocity offsets between features in the small over¬ 
lap region between the FR02 and FR03 data were not misinter¬ 
preted. The small number of double-lined binaries were classi¬ 
fied from inspection of the individual spectra. 

3. Spectral classification 

The spectra were classified by visual comparison to the B-type 
Galactic standards from Sota et al. (2011) and Sana et al. (in 
prep.), while also taking into account the effects of the reduced 
metallicity of the FMC compared to the Galaxy (which can af¬ 
fect the appearance of both the metal and helium absorption 
lines; e.g. Markova et al. 2009). This was achieved for the lu¬ 
minous, low-gravity objects by interpolating between the avail¬ 
able Galactic standards and those for supergiants in the Small 
Magellanic Cloud (Fennon 1997), and with reference to the 
framework developed for the FMC by Fitzpatrick (1988, 1991). 
Classifications of the dwarfs and giants were also informed by 
comparisons with the standards from Walborn & Fitzpatrick 
(1990) and the FMC stars classified by Evans et al. (2006). 

The primary diagnostic lines are the ionisation ratios of sil¬ 
icon, while also taking into account the appearance of the he¬ 
lium and magnesium (A4481) lines. For reference, the primary 
temperature-sequence criteria used to classify the VFTS sample 
are summarised in Tables 1 and 2 (for supergiants and dwarfs, 
respectively). Example temperature sequences for supergiants, 
giants and dwarfs are shown in Figures A.l, A.2 and A.3, re¬ 
spectively. At a given spectral type, luminosity classes were as¬ 
signed from the width of the Balmer lines and (at the earliest 
types) from the intensity of the silicon absorption lines, while 
also taking into account the possible effects of rotational broad¬ 
ening on the spectra. Example luminosity sequences for BOA, 
Bl, and B2.5 types are shown in Figures A.4, A.5, and A.6, re¬ 
spectively. 

The supergiants and bright giants (i.e. class I and II objects) 
have rich absorption-line spectra. Given the importance of spec¬ 
tral peculiarities in the context of CNO abundances, and the in¬ 
fluence of metallicity effects, these luminous objects were classi¬ 
fied independently by three of the authors (CJE, NM, and NRW) 
with a typical consistency of 0.5 subtypes/luminosity class from 
the first pass; final types were adopted following discussion of 
individual objects as required. The larger sample of dwarfs and 
giants, which exhibited a considerable range in data quality, was 
classified by CJE. In many instances these classifications are less 
precise than those possible for the luminous objects, but they 
should provide sufficient morphological information to enable 
robust subdivisions of the sample in future analyses. 

Each merged spectrum was classified in the framework de¬ 
scribed above. After initial classification of all the spectra, care¬ 
ful checks were undertaken to ensure self-consistency within 
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Table 1. Primary temperature-sequence criteria for B-type supergiants. 


Type Criteria 

“BO Si IV AA4089,4116 » Si ill A4553 > He II A4542 

B0.2 Si IV A4089 > Si III A4553 » He II A4542 
B0.5 Si IV A4089 > Si ill A4553; He il A4542 absent/marginal 
B0.7 Si IV A4089 ~ Si III A4553; He II A4686 absent/marginal 
B1 Marginal Si IV A4116; Si ill A4553 > Si IV A4089 
B1.5 Weak/marginal Si IV A4089, weak Mg II A4481 
B2 Si ill A4553 > Mg II A4481; Si IV A4089 absent 

B2.5 Si hi A4553 ~ Mg n A4481; weak Si n AA4128-32 
B3 Si ill A4553 < Mg II A4481; Si II AA4128-32 present 
B5 He IA4471 >MgnA4481; Sill AA4128-32 ~ He IA4121 
B 8 He I A4471 ~ Mg II A4481; Si II A A4128-32 ~ He I A4144 

B9 He IA4471 < Mg IIA4481; Sill AA4128-32 > Hel A4144 

Notes. Example spectra are shown in Figure A.l. 


Table 2. Primary temperature-sequence criteria for B-type dwarfs. 


Type Criteria 

BO Si IV A4089 > Si ill A4553> He II A4542 

B0.2 Si IV A4089 > Si III A4553; He il A4542 marginal 

B0.5 He il A4542 absent; Si ill A4116 marginal/absent 

B0.7 He II A4686 marginal/absent; Si ill A4116 marginal/absent 

B1 Si in A4553 > Si IV A4089; SiivA4116, HeiiA4686 absent 

B1.5 Si IV A4089 absent; Si ill A4553 ~ Mg II A4481 

B2 Si III A4553 < Mg II A4481 

B2.5 Si ill A4553 absent; Si II AA4128-32 marginal 

Notes. Example spectra are shown in Figure A.3. 


the sample. An overlapping sample of ten BO-type stars was 
classified independently (by NRW) for the study by Walborn 
et al. (2014), with excellent agreement (to within one luminosity 
class) in all instances. 

The catalogue of classifications and radial-velocity estimates 
for the B-type objects from the VFTS is presented in Table 7 
(available online). For convenience, supplementary information 
is also provided for each star on its binary status (from Dunstall 
et al., in prep.), Ha morphology (see Section 3.1), optical pho¬ 
tometry (from Paper I), alternative identifications used in the lit¬ 
erature, and previous spectral classifications. 

The S/N ratios of the supergiant spectra were particularly 
good (typically > 100 for both the LR02 and LR03 regions; 
see Figure A.l). This allowed us to provide consistent indi¬ 
cations of the nitrogen enhancement/deficiency compared to 
morphologically-normal stars (via Nstr and Nwk qualifiers), as 
diagnosed from the CNO absorption features in the AA4640- 
4650 region (Walborn 1976). The significant differences in 
CNO morphology are illustrated by the two B1.5 la spectra in 
Figure 1. Atmospheric analysis of the supergiants, including es¬ 
timates of nitrogen abundances and a thorough study of the line¬ 
broadening in their spectra, will be presented by McEvoy et al. 
(2015) and Simon-Diaz et al. (in prep.), respectively. 

Given the wide range in quality of the dwarf and giant spec¬ 
tra, we do not comment on CNO morphologies in their classifi¬ 
cations. Equally, while we employ the ‘n’ and ‘(n)’ indications 
of broadening (see e.g. Table 3 from Sota et al. 2011) at all lumi¬ 
nosity classes, there are likely to be some lacunae (cf. the quan¬ 
titative v e sin % estimates from Dufton et al. 2013). 


3.1. Emission-line objects 

Where appropriate, spectra were initially classified as ‘Be+’ on 
the basis of Fe II emission lines in the LR02 and LR03 spectra 
(see discussion by Evans et al. 2006); the Fe II lines were pre¬ 
ferred over the blue-violet Balmer lines owing to the significant 
(and variable) contamination of the latter by nebular emission. 

Independently of the blue-region classifications, we in¬ 
spected the Ha profiles of every target to investigate their mor¬ 
phologies. The greater resolution of the Ha observations and the 
extended wings from disc emission made it easier to distinguish 
contributions from Be-type Balmer emission and nebular con¬ 
tamination than in the blue, but these classifications should still 
be treated with some caution. The Ha morphologies are encoded 
in the fourth column of Table 7 (with the notation explained in 
the corresponding footnote). 

Objects classified as ‘Be+’ from the LR02/LR03 spectra 
were seen to display Ha emission in all cases. Approximately 
30 additional stars display Be-like Ha emission without Fe II 
emission. The Ha emission in luminous supergiants most likely 
arises from their stellar winds, but for less luminous objects such 
emission is the primary Be diagnostic. To preserve the distinc¬ 
tion of the lines seen in emission, the standard ‘e’ suffix is used 
in the classifications for the (non-supergiant) stars where only 
Ha is seen in emission. 

The complex star-formation history of 30 Dor means that 
our sample will span a range of stellar ages, so we are unable to 
comment on the incidence of the Be-phenomenon for a single, 
co-eval population. Nonetheless, we note that the incidence of 
the Be-type classifications for the dwarf and giant stars is 18% 
(70/388 objects). This is in good agreement with the Galactic 
fraction (> 17%, Zorec & Briot 1997), and that reported for the 
field population near NGC2004 in the LMC (of 16 and 17.5%, 
from Evans et al. 2006 and Martayan et al. 2006, respectively). 

3.1.1. B[e]-type spectra 

Two targets in the Medusa-Giraffe sample show forbidden-line 
emission characteristic of the B[e] category. VFTS 698 is a 
double-lined binary system with a B[e]-like spectrum, compris¬ 
ing what appears to be an early B-type secondary in orbit around 
a veiled, more massive companion (see Dunstall et al. 2012). 
The second B[e] object is VFTS 822, shown in Figure 2 together 
with VFTS 1003 (the B[e]-like object from Paper I). Given the 
intensity of its Si II AA4128-32 and Mg II A4481 absorption, the 
spectrum of VFTS 822 was classified as mid-late B[e] (~B5-8); 
further discussion of its possible pre-main-sequence nature was 
given by Kalari et al. (2014). 

Although the nature of VFTS 1003 remains uncertain at 
present, the number of evolved B[e] stars (at most two) in the 
VFTS is clearly small compared to the number of normal super¬ 
giants. This is in keeping with the small number of sgB[e] stars 
in Galactic clusters (see the appendix of Clark et al. 2013), argu¬ 
ing that the B [e] phenomenon in evolved supergiants must be a 
relatively brief phase or originate from a rare process. 

3.1.2. VFTS 766: A peculiar emission-line star 

VFTS 766 displays strong Ha emission and inspection of its 
blue-region spectrum (see Figure 2) revealed it as a peculiar ob¬ 
ject. At first glance its spectrum appears to be that of a mid-late 
B-type supergiant (B5-B8, as traced by e.g. Si II, Mg II), but 
with superimposed metallic absorption lines, which are consis¬ 
tent with a cooler component. This seemingly composite appear- 
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Wavelength (A) 

Fig. 1. Differing CNO morphologies in early B-type supergiants illustrated by VFTS 431 (‘Nstr’) and VFTS 578 (‘Nwk"). In addition to the 
Balmer lines, the absorption features identified are: C II A4267; He I AA4009, 4026, 4121, 4144, 4388, 4471. 4713; Mg II A4481; N II AA3995, 
4601-07-14-21-31, 4640-43 (blended with O II); O II AA4070, 4254, 4317-19, 4350, 4367, 4414-17. 4591-96, 4650, 4661, 4674-76; Si ill AA4553- 
68-74. 


ance can indicate a circumstellar shell or disc observed edge-on, 
and this hypothesis is supported by the twin-peaked H/3 shell¬ 
like emission profile. Its near- and mid-IR colours 1 place it in 
the region in the J, J— [3.6] colour-magnitude diagram which 
is occupied by early-type stars associated with free-free/bound- 
free emission, and we find no evidence of hot dust (cf. Bonanos 
et al. 2009). 

Similar emission profiles are observed in some low- 
luminosity supergiant B[e] stars (e.g. S59, Gummersbach et al. 
1995), with examples that also appear to lack dusty discs (Gratis 
et al. 2012, although none have yet been discovered in the edge- 
on orientation required to yield shell profiles). However, we 
speculate that VFTS 766 is a classical Be star, in which the twin- 
peaked emission in the wings of the Balmer profiles gives rise to 
the supergiant-like appearance. An example of this effect is the 
development of a circumstellar disc in the classical edge-on Be 
star o And, which conspires to yield a narrower absorption fea¬ 
ture than the underlying photospheric profile (Clark et al. 2003). 
Another example is provided by the Be star 48 Lib (see Rivinius 
et al. 2013, and references therein), which has previously been 
classified as a supergiant (Houk & Smith-Moore 1988). 

This possible explanation for the nature of VFTS 766 is sup¬ 
ported by its magnitude relative to other stars in our sample. 
Spectroscopic monitoring to identify an episode of disc loss, dur¬ 
ing which the (uncontaminated) photospheric Balmer line pro¬ 
files might be observed, would clarify its luminosity type and 
true evolutionary status. 

3.1.3. Blue hypergiants 

We identifed five blue hypergiants (BHGs) in our sample, 
namely: VFTS 003 (B1 Ia+), 424 (B9 Ia+p), 430 (B0.5 
Ia+((n)) Nwk), 458 (B5 Ia+p) and 533 (B1.5 Ia+p Nwk). Recent 
theoretical and observational studies suggest that early-type 
BHGs might be massive stars which have passed through a blue 

1 VFTS 766: J= 15.14, 77=15.01, K s = 14.76 (from Paper I); 
[3.6] = 14.28, [4.5] = 14.06. [5.8] = 13.64 (from Meixner et al. 2006). 


supergiant phase and are about to enter a luminous blue variable 
(LBV) phase (Clark et al. 2012; Groh et al. 2014). This suggests 
that VFTS 003, 430, and 533 might be on the verge of LBV- 
like behaviour, although we note that VFTS 430 appears some¬ 
what subluminous compared to the blue supergiants in 30 Dor, 
and is significantly fainter than VFTS 533 (even accounting for 
its apparently greater line-of-sight extinction, see photometry 
in Table 7). Later-type BHGs have been observed to undergo 
LBV-like photometric and spectroscopic behaviour (e.g. Smith 
et al. 2004; Clark et al. 2012, with potential implications for 
VFTS 424 and 458 here), suggesting that these two evolution¬ 
ary phases are essentially synonymous for such stars. 

On further inspection of their spectra, broad emission wings 
were identified in the H/3 profiles of each of the BHGs, except for 
VFTS 430. This Hi morphology is seen in some LB Vs and other 
luminous supergiants (e.g. Walborn & Fitzpatrick 2000), and is 
also present in VFTS 739 in our sample (classified as A0 Ip but, 
given its omission from Paper I, included here in the B-type 
sample). The spectra and evolutionary status of these potential 
LBV precursors/candidates will be discussed in more depth by 
Walborn et al. (in prep.). 


3.1.4. Comparisons with X-Ray observations 


We compared our sample with the point-source detections from 
observations in 30 Dor with the Chandra X-ray Observatory 
(Townsley et al. 2006, 2014); none of the VFTS B-type stars 
is detected in the extant X-ray data. A new X-ray Visionary 
Project is now underway, ‘The Tarantula Revealed by X-rays (T- 
ReX)’, which will observe 30 Dor with Chandra for 2 Ms. This 
will push far deeper than the existing data, and we will revisit 
the VFTS sample once the T-ReX observations and analysis are 
complete. 
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Wavelength (A) 

Fig. 2. Spectra of three peculiar Be-type stars from the VFTS, in which each spectrum has been smoothed and rebinned to R = 4 000 for clarity. 
In addition to the Balmer lines, the emission lines identified in VFTS 1003 are [S n] A4069; Fe II AA4173. 4179, 4233, 4303, 4352, 4385, 4481, 
4491, 4508. 4515, 4520-23, 4549, 4556; [Fe II] AA4244, 4277, 4287, 4358-59, 4414-16. Emission lines identified in VFTS 822 are Fe II AA4584, 
4620.5, 4629. Absorption lines identified in VFTS 766 are He I AA4009, 4026, 4121, 4144, 4388, 4471, 4713; Si II AA4128-32; Mg II A4481. 


4. Stellar radial velocities 

4.1. Single stars 

Radial-velocity estimates for the stars identified as apparently 
single by Dunstall et al. (in prep.) were obtained from Gaussian 
fits to selected absorption lines in the merged LR02 and LR03 
spectra of each object (at the native resolution of the data, with¬ 
out the degrading applied for classification purposes). This ap¬ 
proach is consistent with that used by Sana et al. (2013) for 
analysis of the O-type stars, albeit (necessarily) using a different 
subset of absorption lines. We note that the choice of Gaussian 
profiles may not be appropriate for fitting the wings of the lines 
in the presence of different broadening mechanisms, but they 
should suffice for robust estimates of the centres of our chosen 
lines. 

We also considered a cross-correlation method to estimate 
the radial velocities. However, given the range in temperature 
and luminosity of the sample, as well as the varying data quality 
and problems of nebular contamination, the selection of suitable 
template spectra for cross-correlation was not trivial and would 
have introduced additional uncertainties (relating to e.g. the pro¬ 
jected rotational velocities, v e sin i). Thus, we did not pursue this 
approach further. 

4.1.1. Line selection 

The large range in v e sin i and effective temperature for the sam¬ 
ple meant that it was not possible to use a single set of diagnostic 
lines for all stars. Moreover, the large intrinsic width and nebu¬ 
lar contamination of the Balmer lines precluded their use; sim¬ 
ilarly, nebular emission affects some of the stronger He I lines. 


Following the approach taken by Dufton et al. (2013), we used 
different sets of diagnostic lines (as summarised in Table 3) de¬ 
pending on the v e sin i of the star. 

For narrow-lined spectra (defined as v e sin i < 150 km s -1 ), 
we used three non-diffuse He I lines and four isolated metal lines 
listed as ‘Set 1’ in Table 3. Where possible, radial velocities were 
estimated from Si III A4553 in both the merged LR02 and LR03 
spectra, giving up to eight individual estimates for the narrow- 
lined objects. 

We note that He I A4121 is potentially affected by blending 
with O II AA4120.28/54. In many instances these features could 
be resolved and a double Gaussian profile was fitted to estimate 
the He I and O II line centres; in cases where the O II lines were 
weaker or unresolved, they do not appear to have unduly influ¬ 
enced the results (see Table 3, in which the residual for He I 
A4121 is consistent with the other lines). 

Radial-velocity estimates for stars with broader lines were 
more difficult, particularly for spectra with significant nebular 
contamination. In these cases, velocities were estimated by fit¬ 
ting the five He I lines listed as ‘Set 2’ in Table 3, for which the 
profiles are expected to be effectively symmetric 2 . 

Representative examples are shown in Figure 3. The lower 
panel shows fits to the Si III triplet for two narrow-lined stars, 
VFTS 053 and 159, which have v e sin * < 40 (i.e. less than the 
velocity resolution of the data) and 96 ±9 km s -1 , respectively 
(from Dufton et al. 2013). As an example of a spectrum of a 
rapidly-rotating star which also suffers nebular contamination. 


2 He I 3 P —> 3 D transitions were excluded because of the presence of 
significant ! P —> 3 F components. 
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Table 3. Rest wavelengths used to estimate radial velocities (v r ) for the 
narrow- and broad-lined stars (Sets 1 and 2, respectively). Values in the 
final column are average differences (A) between the velocity estimates 
for each line and the final velocity for each star. 


Ion 

A [A] 

A [km s -1 ] 

He I 

4120.82 

Set 1 

-0.8 ±6.7 

CII 

4267.13 

1.1 ±5.1 

He I 

4437.55 

-0.4 ±5.4 

Si III 

4552.62 

0.3 ±6.1 (LR02) 

Si III 

4567.84 

— 1.6 ±5.8 (LR03) 
-0.8 ±6.4 

Si III 

4574.76 

— 1.7 ± 7.0 

He I 

4713.15 

2.1 ±4.5 

He I 

4009.26 

Set 2 

—0.8 ± 11.3 

He I 

4143.76 

0.1 ±9.0 

He I 

4387.93 

2.7 ±9.1 

He I 

4713.15 

2.8 ±7.1 

He I 

4921.93 

-3.5 ±8.7 


Notes. Wavelengths are taken from the NIST atomic spectra database 
(Kramida et al. 2012); the adopted value for C II is the average of the 
two stronger transitions. 


the upper panel of Figure 3 shows the Gaussian fit to He I A4388 
for VFTS 636 (v e sin * = 371 ± 35 km s" 1 ; Dufton et al. 2013). 

For three stars with later/peculiar types - VFTS 272 (possi¬ 
ble shell star), 739 (AO Ip), and 766 (B5-8e, see Section 3.1.2) - 
we employed weak Fe II and Ti II absorption lines to estimate ve¬ 
locities, although we note that in VFTS 272 and 766 they might 
not be representative of the true (stellar) radial velocity, as these 
lines could trace cooler features associated with a disc of mate¬ 
rial (see e.g. the discussion of VFTS 698 by Dunstall et al. 2012). 




Wavelength (A) 

Fig. 3. Examples of Gaussian fits used to estimate radial velocities (ti r ) 
for the B-type sample. Upper panel: Fit to He I A4388 in a rapidly- 
rotating object (VFTS 636) with nebular contamination. Lower panel: 
Fits to the Si ill triplets (AA4553-68-74) of VFTS 053 and 159. 


4.1.2. Mean velocities 


Radial-velocity estimates for the single stars were obtained by 
calculating an average of the n individual measurements (v,). 
weighted by the estimated central depth (di) from the Gaussian 
fit of each line, i.e. 


v r 


EE 1 v i d i 

EEl^ 


( 1 ) 


The sample standard deviation (a) for each star was similarly 
calculated as a weighted average: 


2 = = 

EEl di 


( 2 ) 


In cases where the central depths of the absorption lines were 
less than twice the inverse of the continuum S/N ratio, the in¬ 
dividual estimates were excluded from the above calculations 
(to avoid their larger uncertainties degrading the results from 
stronger lines). For stars with estimates from three (or more) 
absorption lines, we present their estimated mean velocities and 
standard deviations in columns 7 and 8 of Table 7. For complete¬ 
ness, the line-by-line estimates from Sets 1 and 2 for the single 
stars are detailed in Tables A. 1 and A.2, respectively; the mean 
uncertainty of the estimates for the 307 (presumed) single stars 
is 7.8 km s^ 1 (and with a median uncertainty of 6.2 km s^ 1 ). 

As a check on our results, we used the methods described 
above to analyse a sample of the late O-type stars from the 


VFTS, finding excellent agreement between our velocity es¬ 
timates and those from Sana et al. (2013, see discussion in 
Section A.2.1). We also investigated the internal consistency of 
our results from the selected lines by calculating the average dif¬ 
ference (A) of the individual measurements for each line com¬ 
pared to the adopted mean velocity for each target. The means 
and standard deviations of these differences are given in col¬ 
umn 3 of Table 3, with no significant shifts seen for any of the 
diagnostic lines compared to the adopted radial velocities. 

4.2. Radial velocities of the binaries 

To complement the analysis of the single stars, we also inves¬ 
tigated the radial velocities of the single-lined binaries in the 
sample. We obtained single-epoch velocity estimates (tingle) for 
each system using measurements from the best LR02 spectra (in 
terms of S/N ratio of the co-added exposures from a given night, 
which span three hours at most). 

Given the reduced spectral coverage of the LR02 data com¬ 
pared to the full LR02 + LR03 range, we used all of the available 
lines from Table 3 to estimate tangle (except for He I A4121, 
given the nearby O II feature and the lower S/N ratio of the 
single-epoch spectra compared to the quality of the co-added 
spectra of the single stars). Line-by-line estimates and the date 
for which they were measured are listed in Table A.3. Mean ve- 
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locities and standard deviations for the selected epoch were then 
estimated using the same methods as for the single stars. 

The mean single-epoch velocities serve as an absolute ref¬ 
erence point for the inter-epoch differential velocities calculated 
by Dunstall et al. (in prep.) yielding absolute estimates for each 
epoch. We then calculated the mean and standard deviation of 
the multi-epoch estimates to obtain a velocity estimate for each 
single-lined binary. These values are presented in Table 7, and 
are flagged as ‘B’ (binary) in the ninth column. We note that the 
typical uncertainty on the cross-correlation results from Dunstall 
et al. was less than 5 km s' 1 , while the typical error on the ab¬ 
solute estimates determined here from the single-epoch data was 
9 km s' 1 ; i.e. the uncertainties in our estimates will generally be 
dominated by the analysis presented here. 

Full orbital solutions are presently unavailable for the B-type 
binaries, so we strongly caution the reader that the velocity esti¬ 
mates in Table 7 are not their centre-of-mass, systemic velocities 
(although, as discussed in the Appendix, in many instances they 
may provide reasonable estimates). 

4.3. Spatial distribution 

In addition to NGC2070 (the main 30 Dor cluster, which in¬ 
cludes R136 at its core) and NGC2060 (the association to the 
south-west), there are two smaller clusters in the VFTS sur¬ 
vey: Hodge301 (Hodge 1988) and SL693 (Shapley & Lindsay 
1963); their relative locations are shown in Figure 4. Both clus¬ 
ters appear older than NGC2060 and 2070 as they contain 
evolved supergiants, from B- through to M-types. 

We calculated weighted mean velocities (fy. and their stan¬ 
dard deviations) for the full sample of single stars, for the 
subsamples associated with the four distinct stellar assocations 
in the survey region (NGC2060, NGC2070, Hodge 301, and 
SL639, as defined in Table 4), and for the field population of 
stars outside the four clusters. These results, which exclude can¬ 
didate runaway stars (see Section 6), are summarised in Table 4. 

The results for the two younger star-forming regions (i.e. 
NGC 2060 and 2070) are in excellent agreement with the field 
population of B-type stars in the region, but the two older clus¬ 
ters (Hodge 301 and SL639) appear at slightly lower recession 
velocities, suggesting that they are kinematically distinct from 
the rest of the sample. The weighted mean velocity for the 
273 single stars excluding these two clusters (and potential run¬ 
aways) is uf = 271.6 ± 12.2 km s' 1 . 

5. Spectral content of the older clusters 

The VFTS data comprise the first comprehensive spectroscopy 
in Hodge 301 and SL639, so we briefly discuss these clusters in 
turn. 

5.1. Hodge 301 

There are 20 VFTS targets within a 20" radius (chosen to conser¬ 
vatively identify likely members, and corresponding to ~4.9 pc, 
assuming a distance modulus of 18.5 mag). Fifteen of these were 
classified as B-type and their likely membership of the cluster is 
indicated in the final column of Table 7 3 . 


3 The five others are: VFTS 310 (09.7 V:, Walborn et al. 2014) and 
VFTS 271, 281, 289, and 294 (A7 II, Mid-late K, Late G/Early K, and 
A0 lb, respectively, from Paper I). Three additional stars (VFTS 263, 
291, and 317) are also nearby, at radii from the adopted centre of ~25". 


From analysis of Hubble Space Telescope (HST) imaging, 
Grebel & Chu (2000) estimated Hodge 301 to be 20-25 Myr 
old. From narrow-band imaging they also identified 19 Be-type 
stars, eight of which (VFTS 272, 276, 279, 282, 283, 287, 293, 
and 301) have spectroscopy from the VFTS. The Ha spectra of 
each of these contain significant nebular emission (with mul¬ 
tiple components in all sightlines except towards VFTS 293). 
Nonetheless, broad Be-like Ha emission is present in each ob¬ 
ject, except for VFTS 276 (GC00-Be2). 

Grebel & Chu (2000) noted two potential blue stragglers in 
the cluster, which they argued were too luminous for the main- 
sequence population (but less luminous than expected for su¬ 
pergiants). The first, VFTS 270 (WB97-2), is a seemingly unre¬ 
markable B3 lb star. The second is VFTS 293 (WB97-9, GC00- 
Bel), classified here as B2 III-II(n)e. From two previous ob¬ 
servations, Walborn & Blades (1997) noted small shifts of its 
absorption features (relative to superimposed nebular emission) 
and possible weak He II A4686 emission at one epoch. If the 
He II emission were real, Walborn & Blades speculated that it 
might be an X-ray binary. There is no evidence of A4686 emis¬ 
sion in the LR03 FLAMES spectra (albeit all taken on the same 
night, see Paper I), but small line-profile variations are seen in 
the absorption features in the multi-epoch LR02 spectra. 

More intriguing is VFTS 310 on the periphery of the clus¬ 
ter and classified as 09.7 V: (Walborn et al. 2014). With such a 
classification, VFTS 310 would be expected to be significantly 
younger than the other cluster members, as shown by its loca¬ 
tion in the Hertzsprung-Russell (H-R) diagram of the cluster 
(Figure 5, see Section 5.3 for details). The star appears oth¬ 
erwise unremarkable, with a low projected rotational velocity 
(v e sin? < 40km s' 1 , Dufton et al. 2013) and a radial veloc¬ 
ity (v r = 272.3 ±2.8 km s' 1 , Sana et al. 2013) which is con¬ 
sistent with those in NGC 2070 and the local field population. 
Moreover, it does not appear to be a runaway from Hodge 301 
(see Section 6). Whether this is a tme blue straggler or simply a 
line-of-sight coincidence remains uncertain at this time. 

5.2. SL 639 

Adopting a 20" radius (~4.9 pc), sixteen of the VFTS targets 
are within SL693. Fourteen are in the B-type sample, with their 
membership of SL 639 indicated in the final column of Table 7 4 . 
Ten of these stars are early-B dwarfs/giants, with six of them dis¬ 
playing Be-like Ha emission. The remaining four are comprised 
of two supergiants (VFTS 827, B1.5 lb; VFTS 831, B5 la) and 
two bright giants (VFTS 826, B1 Iln; VFTS 829, B 1.5-2 II). 

There are previous age estimates for SL693. From consider¬ 
ation of its main-sequence turn-off, Hodge (1983) estimated its 
age to be 18 Myr (although he noted a large uncertainty com¬ 
pared to the theoretical models of the time). A larger estimate 
of 30 Myr was given by Santos et al. (1995) from its reddening- 
corrected, integrated (U — B)q colour and adopting an empirical 
calibration of colour vs. age from Bica et al. (1990). 

5.3. H-R diagram 

Motivated by the apparent similarities of the two clusters in 
terms of their spectral content, we constructed an H-R dia- 

4 The two others are VFTS 820 and 828 (A0 la and Early M, respec¬ 
tively. from Paper I). The only other VFTS target within a 35" (8.5 pc) 
radius is VFTS 839 at ~25" (6pc), classified as G-type in Paper I. All 
three of these cooler objects have radial velocities consistent with mem¬ 
bership of the LMC. 
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RA (J2000) 


Fig. 4. Spatial distribution of the B-type single stars (cyan circles) and spectroscopic binaries (open blue squares) with radial-velocity estimates; 
the nine candidate runaways, labelled by their VFTS identifiers, are marked in red. The spatial extents of NGC 2070, NGC 2060. SL 639, and 
Hodge 301 (as defined in Table 4), are indicated by the overlaid dashed circles. The underlying image is from a V-band mosaic taken with ESO’s 
Wide Field Imager on the 2.2 m telescope at La Silla. 


gram in order to compare the properties of their likely members. 
Ahead of detailed quantitative analysis of the VFTS spectra and 
pending extinction analysis of the whole sample, we adopted ef¬ 
fective temperatures (T e g) on the basis of the spectral types in 
Table 3 using the calibrations of Trundle et al. (2007), and in¬ 
terpolating between the Galactic and SMC temperatures for A- 
type supergiants from Evans & Howarth (2003); cooler objects 
were omitted given the larger uncertainties in their spectral clas¬ 
sifications. Stellar luminosities were estimated using the opti¬ 
cal photometry from Paper I, intrinsic colours from Fitzgerald 
(1970), bolometric corrections from Balona (1994), and assum¬ 
ing Ry = 3.5 (see Appendix C of Doran et al. 2013). 


The H-R diagram is shown in Figure 5, with SL639 and 
Hodge 301 members indicated by the filled black and open red 
circles, respectively. Evolutionary tracks and isochrones (for 
non-rotating models) from Brott et al. (2011) are also plotted for 
comparison. As noted earlier, the late O-type star in Hodge 301, 
VFTS 310, stands out as a blue straggler with T e g > 30 kK. 

Tailored atmospheric and extinction analyses of these ob¬ 
jects should help to provide further insight in due course. 
Nonetheless, the stars in SF639 appear marginally younger 
than those in Hodge 301. From qualitative comparison with the 
isochrones, we estimate an age in the range of 10-15 Myr for 
SF639, significantly younger than the estimate of 30 Myr from 
Santos et al. (1995). There is a larger spread in the inferred ages 
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Table 4. Mean radial velocities ( v r ) and standard deviations (cr) of the stellar associations identified within 30 Doradus. 


Region 

No. of stars 

Radius 

['] 

a 8 

[J2000.0] 

Vr ± (7 

All 

298 




270.4 ± 12.4 

NGC 2070 

78 

2.40 

05 3842.0 

-690602.0 

271.2± 12.4 

NGC 2060 

22 

2.40 

05 3745.0 

-691015.0 

273.7 ± 13.4 

Hodge 301 

14 

0.33 

05 3817.0 

-690401.0 

261.8± 5.5 

SL 639 

11 

0.33 

05 39 39.4 

-691152.1 

253.5 ± 3.9 

Field 

173 




271.3± 11.8 

All (excl. H301/SL639) 

273 




271.6± 12.2 


Notes. These statistics exclude the nine candidate runaways identified in Table 5. For consistency, the central positions adopted for NGC 2060 and 
NGC2070 are those used by Sana et al. (2013). Approximate centres for Hodge 301 and SL639 were adopted using the coordinates of WB97-5 
(Walborn & Blades 1997) and VFTS 828, respectively. 



Fig. 5. Hertzsprung-Russell diagram for the VFTS objects observed in 
SL639 (black circles) and Hodge 301 (red, open circles). Also plotted 
are evolutionary tracks for initial masses of 8, 10, 12, 15, 20, 25 Mq 
( solid lines), the zero-age main sequence (dashed line), and isochrones 
(3, 5, 10, 15, 20, and 25Myr, dotted lines) for the non-rotating mod¬ 
els from Brott et al. (2011). The apparent blue straggler in Hodge 301 
(VFTS 310) is a late O-type star. 


for the stars in Hodge 301 of 15 ± 5 Myr, overlapping with the 
lower estimate of 20 Myr from the photometric study of Grebel 
& Chu (2000). Adopting isochrones calculated for the rotating 
models from Brott et al. had little impact on these results; e.g. 
assuming an intial rotational velocity of 301 km s , the inferred 
ages were only marginally younger (by ~1 Myr). 

6. Candidate runaways 

6.1. Single stars 

To identify candidate runaways in the B-type sample we com¬ 
pared the radial-velocity estimates for each star with the initial 
weighted mean and standard deviation calculated for the stars in 


NGC 2060, NGC 2070, and the surrounding field population (i.e. 
excluding Hodge 301 and SL639). Employing a threshold of 
Sv r > 3cr we identified eight candidate runaways and then re¬ 
calculated the mean velocities (and dispersions) excluding these 
objects. Reapplying a 3 -cr threshold provided one additional can¬ 
didate, while a second iteration did not reduce the velocity dis¬ 
persion further. The mean velocities and dispersions for the dif¬ 
ferent spatial samples (excluding potential runaways) are pre¬ 
sented in Table 4. The cumulative distribution of the velocity 
estimates for the 307 single stars is shown in Figure 6, with the 
±3cr threshold (i.e. ±36.6 km s _1 ) for identification of runaway 
stars indicated by the vertical dotted lines. 

Our adopted velocity threshold is comparable in magnitude 
to the limit used by Blaauw (1961), which was that candidate 
runaways have a space velocity which differs by more than 
40 km s -1 from the systemic velocity of the region in which they 
formed. (In this context, our radial-velocity estimates should 
serve as lower limits to the 3-D space velocities.) In contrast, 
Portegies Zwart (2000) adopted a lower threshold of 25 km s -1 
in his theoretical study of the origins and characteristics of run¬ 
aways. Given the limitations of the B-type data and the relatively 
small number of lines available for some objects, we did not con¬ 
sider a comparably low threshold. 

Details of the nine candidate runaways are summarised in 
Table 5, and their locations are shown in Figure 4. Four of 
the candidates have differential velocities in the range of 3.2- 
3.6 cr, thus their status as candidates will depend strongly on the 
adopted velocity criteria. Moreover, our adopted velocity thresh¬ 
old (the dispersion of the weighted mean velocities in the sam¬ 
ple) was employed without consideration of the uncertainty on 
the radial-velocity estimate for a given star. Thus, VFTS 467, al¬ 
though formally with a radial velocity of 4.2<r from the systemic 
value, has a large uncertainty on its estimated velocity. Indeed, 
only three He I lines were available for this star and the A4144 
estimate (277.5 km s _1 ) is comparable to the systemic value. 

The choice to use weighted means to investigate the 
spatially-different samples in Section 4.3 was motivated by the 
fact that approximately 10% of our velocity estimates have un¬ 
certainties of > 15 km 1 (each of which is a broad-lined star, 
i.e. used the Set 2 lines from Table 3). To investigate the de¬ 
pendence of the number and identity of the candidate runaways 
on these adopted means, we repeated the above steps for both 
the (unweighted) mean and median velocities of the sample. 
The mean velocity for the 281 stars (not including those in 
Hodge 301 and SF639) was 269.8 = 16.4km s -1 , with a median 
value of 270.7 km s _ . Adopting either of these systemic values 
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Table 5. Candidate runaway stars identified from the radial-velocity analysis of the VFTS B-type spectra. 


VFTS 

Classification 

Vi 

[km s' 1 ] 

(v r - v r ) 

[km s' 1 ] 

| V r — Vr\!a 

v e sin i 

[km s' 1 ] 

['] 

Vd 

[pc] 

242 

B0 V-IV 

216.0± 3.6 

-55.6± 12.7 

4.6 

<40 

5.96 

86.6 

298 

Bl-2 V-IIIe± 

332.1 ± 4.0 

60.5 ±12.8 

5.0 

431 ± 15 

3.92 

57.1 

358 

B0.5: V 

164.7 ± 10.7 

—106.9 ± 16.2 

8.8 

345 ± 22 

3.02 

44.0 

368 

Bl-3 V 

229.6 ± 12.4 

-42.0 ± 17.4 

3.4 

380: 

1.20 

17.5 

467 

Bl-2 Ve± 

322.3 ±30.2 

50.7 ±32.6 

4.2 

355 ± 29 

9.54 

138.7 

616 

B0.5: V 

227.4 ± 10.7 

-44.2 ± 16.2 

3.6 

<40 

1.26 

18.3 

625 

B1.5 V 

310.7 ± 9.2 

39.1 ±15.3 

3.2 

64 ± 12 

2.57 

37.4 

831 

B5 la 

210.4 ± 6.2 

61.2 ± 13.7 

5.0 

41 

7.91 

115.0 

851 

B2III 

229.4 ± 0.3 

-42.2 ± 12.2 

3.5 

<40 

7.18 

104.4 


Notes. Estimates of v e sin* are from Dufton et al. (2013) except for VFTS 831. which is taken from McEvoy et al. (2015). Projected distances 
from the core of R136 are from Paper I. 


and a 3-d threshold (taking the dispersion as 16.4 km s' 1 ) led to 
identification of the five most significant runaways from Table 5 
(i.e. VFTS 242, 298, 358, 467, and 831), with a second iteration 
providing one additional candidate (VFTS 616). The final mean 
velocity from these calculations is 270.1 ± 13.9, in good agree¬ 
ment with the equivalent subsample in the final row of Table 4. 
From these checks we conclude that the four most significant 
candidates (VFTS 242, 298, 358, and 831) are secure as poten¬ 
tial runaways, with five more tenative candidates as discussed 
above. 

For the members of Hodge 301, VFTS 309 is a poten¬ 
tial outlier (at 3.6d) compared to the mean velocity for 
the cluster in Table 4. Indeed, excluding its i> r estimate 
from the calculated mean and dispersion for Hodge 301 gives 
*V =261.5 ± 5.0 km s' 1 , such that VFTS 309 is then a 4.1 cr out¬ 
lier. However, its velocity is also in reasonable agreement with 
those for stars in NGC 2060 and NGC 2070, so we are unable to 
confirm it as a potential runaway from Hodge 301. Finally, we 
note that the velocity estimate for VFTS 310 is not significant as 
a potential runaway from Hodge 301. 

With our adopted velocity criteria, the nine candidate run¬ 
aways correspond to 2.9% of the sample of 307 single stars 
with radial-velocity estimates. This lower fractional share of run¬ 
aways compared to the O-type stars (Sana et al., in prep.) is 
in qualitative agreement with recent N-body predictions for the 
mass-spectrum of runaways (Banerjee et al. 2012; Banerjee & 
Kroupa 2012), although quantitative comparison with the pre¬ 
dictions is complicated by the different criteria employed (e.g. 
Banerjee et al. used a radial distance threshold for a star to be 
considered a runaway). We note that the distribution in Figure 6 
is slightly asymmetric, with two thirds of the candidates hav¬ 
ing negative differential velocities. This could point to veiling of 
more distant runaways by the nebula, although this is not seen 
in the results for the O-type stars (in which a greater fraction 
are seen to have positive differential velocities, Sana et al., in 
prep.), suggesting that small-number statistics likely explain the 
asymmetry. 

6.2. Rotational velocities of candidate runaways 

Estimated projected rotational velocities (from Dufton et al. 
2013; McEvoy et al. 2015) for our candidate B-type runaways 
are included in Table 5. There appears to be a dichotomy in their 
rotation rates, as highlighted by Fig. 7. Five stars have small ro¬ 
tational velocities, with the estimates for four stars effectively 
limited by the velocity resolution of the data (i.e. ^40 km s' 1 ) 



175 200 225 250 275 300 325 

v T (km s 1 ) 


Fig. 6. Cumulative distribution of estimated radial velocities for the 
307 single B-type stars. The dashed vertical line indicates the weighted 
mean velocity of the sample. v r =271.6km s' 1 (calculated excluding 
stars in the two older clusters and the candidate runaways). The vertical 
dotted lines are the 3 a thresholds adopted for identification of potential 
runways (±36.6 km s' 1 ). 


and the fifth with v e sin* = 64± 12km s' 1 . Meanwhile, the re¬ 
maining four stars each have estimates of v e sin i > 345 km s' 1 . 
To investigate the significance of these results, we ran Monte 
Carlo simulations using the v e sin i distributions from Dufton 
et al. (2013). We note that the results from Dufton et al. were for 
the unevolved B-type sample (i.e. dwarf and giant objects), so 
we exclude VFTS 831 from the following discussion. 

We generated random rotational velocities for eight stars 
drawn from the distribution of Dufton et al., and repeated 
this 10 6 times to estimate the probability of recovering rota¬ 
tional velocities consistent with our eight (non-supergiant) run¬ 
away candidates (i.e. four stars with v e sin* < 65, and four with 
v e sin* > 345 km s' ). We then repeated these runs 10 4 times 
to estimate the dispersion on the probabilities. In a second test 
we limited the low velocity subsample to the three stars with 
*; e sin* < 40km s' 1 (i.e. excluding VFTS 625, which has the 
least significant runaway velocity), together with the four rapid 
rotators. We found low probabilities for recovering the observed 
v e sin * values from both tests (as summarised in Table 6), with 
significance levels of 3.0 and 3.2d, respectively. 
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Table 6. Probability tests to investigate the rotational velocities of the 
candidate runaway stars compared to those for the full (unevolved) B- 
type sample from Dufton et al. (2013). 


Test (from 10 4 repeats of 10 b runs) 

n 

P 

4x» e sin i< 65 & 4 x u e sin i > 345 km s _1 

10 ±3 

0.0010 

3 x v e sin * < 40 & 4 x u e sin i > 345 km s _1 

6±2 

0.0006 

4 x v e < 100 & 4 x v e > 340km s~ 4 

395 ± 20 

0.0395 

3x» e < 60 & 4 x v e > 340km s~ 4 

181 ±13 

0.0181 


Notes, ‘n’: number of Monte Carlo runs (from 10 6 sets) for which the 
simulated rotational velocities are consistent with those of the candidate 
runaways (and translated into a fractional probability, P, in the final 
column); the quoted uncertainties are the standard deviations from 10 4 
repeat runs. 


500 


400 

m 

a 300 
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m 
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> 
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0 20 40 60 80 100 120 

<5v r | (km s _1 ) 

Fig. 7. Comparison of the projected rotational velocities ( v e sin i) with 
the differential radial velocities (<5w r ) compared to the mean systemic 
velocity of the B-type sample. Uncertainties on v e sini of ±10% are 
adopted for VFTS 368 and 831 (with the latter plotted as an open sym¬ 
bol given its supergiant classification). 





We also investigated the probabilities compared to the cumu¬ 
lative distribution function of unprojected rotational velocities 
(v e ) from Dufton et al. (their Table 6). In this case we are limited 
by the unknown inclination angles of our targets but, assuming 
some relatively conservative lower-limits, we repeated tests sim¬ 
ilar to those above (in which we adopt the nearest available value 
from Dufton et al. for our upper limit, i.e. 340km s~ 1 ). These 
tests are formally less significant (Table 6), but provide further 
support that the rotational velocities of our candidate runaways 
do not appear to be randomly drawn from the distribution for 
the larger B-type sample. We note that this agrees with similar 
conclusions regarding the rotational velocities of the O-type run¬ 
aways (Walborn et al. 2014, Sana et al., in prep.). 

6.3. Atmospheric analysis of two candidate runaways 

Atmospheric analysis of the B-type spectra is now underway 
within the VFTS project. In the context of the current study, we 
were particularly motivated to investigate the properties of the 
two most noteworthy runaway candidates, VFTS 358 and 831; 
we now comment briefly on our findings for these two objects. 


6.3.1. VFTS358 

This star has the largest differential velocity from the 
estimated systemic value (Sv T — —106.9 ± 16.2km s _1 ) and 
also has a relatively large projected rotational velocity 
( v e sin * = 345 ± 22 km s _1 ; Dufton et al. 2013). This prompted 
us to investigate its physical properties, in part to see if there was 
any evidence of past binary interaction (see e.g. de Mink et al. 
2014). 

Detailed analysis was complicated by the rapid rotation, 
which made estimation of the atmospheric parameters and el¬ 
ement abundances difficult. However, using similar methods as 
those used by Hunter et al. (2007), we undertook a preliminary 
analysis using a grid of model atmospheres calculated using the 
TLUSTY and SYNSPEC codes (Hubeny 1988; Hubeny & Lanz 
1995; Hubeny et al. 1998; Lanz & Hubeny 2007); further details 
of the grids are given by Ryans et al. (2003) and Dufton et al. 
(2005) 5 . 

We estimated the effective temperature (7(, ff ) and gravity 
(log g) of VFTS 358 by fitting rotationally-broadened profiles to 
the H(5, H 7 and He II A4686 lines, assuming solar helium abun¬ 
dances. This gave estimates of T e g ~ 29 kK and log g ~ 3.5 dex, 
although we were unable to estimate the microturbulence be¬ 
cause of the paucity of observable metal lines (which also pre¬ 
cluded use of the silicon ionisation balance to estimate the tem¬ 
perature). 

Absorption from N II A3995 appears to be present in the 
spectrum (with an equivalent width of ~ 120 mA), although its 
identification is marginal, as shown in Figure 8 . Rotationally- 
broadened model spectra for our adopted physical parameters 
(and an assumed microturbulence of 5 km s _1 ) are shown for 
three nitrogen abundances: 12 + log(N/H) = 6.9 (the approxi¬ 
mate baseline N-abundance for the LMC, see discussion by 
Hunter et al. 2007), 7.7, and 8.5 dex plotted in red, green, and 
blue, respectively. This comparison suggests significant nitro¬ 
gen enrichment, with an abundance of > 8.5 dex (the maximum 
abundance considered in our grid). Tests for theoretical pro¬ 
files with a microturbulence of 20 km s" 1 (which is larger than 
normally found for stars with this gravity, e.g. Hunter et al. 
2007) implied a similarly large nitrogen abundance. Fitting a 
rotationally-broadened profile to the N II A3995 line leads to 
v e sin i = 342 km s -1 (cf. 345 km s -1 from the helium lines) and 
a central wavelength of 3995.06 A (cf. the laboratory wavelength 
of 3995.00 A). This agreement provides some support that the 
feature is indeed from N II absorption. 

The neutral helium lines in the observed spectrum appear 
somewhat stronger than those in the best-fitting model (e.g. the 
He I AA4009, 4026 lines in Figure 8 ) suggesting that the he¬ 
lium abundance might also be enhanced. In turn, this would 
tend to decrease the T e g estimated from the He II line (and the 
surface gravity). Tests using cooler models indicated that the 
change in nitrogen abundance would be modest. For example, 
a change in effective temperature of — 2 kK leads to a change in 
12 + log(N/H) of approximately —0.2 dex. Hence, the nitrogen 
abundance would still be significantly greater than the baseline 
value for the LMC. 

At present we are limited by the quality of the VFTS spec¬ 
trum but from this preliminary analysis there is evidence of sig¬ 
nificant nitrogen enrichment (and maybe also helium) in the pho¬ 
tosphere of VFTS 358. Indeed, the estimated nitrogen abundance 
is the largest obtained so far from the VFTS survey (including 
results from the B-type supergiants from McEvoy et al. (2015). 

5 See also http://star.pst.qub.ac.uk 
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Wavelength (A) 


Fig. 8. Section of the VFTS 358 spectrum showing, from left-to-right, 
the He, N II A3995, and He I AA4009, 4026 lines; we note that He is 
complicated by absorption from the interstellar Ca H line (as well as 
nebular emission). Model spectra for our adopted physical parameters 
and 12 + log (N/H) = 6.9, 7.7, and 8.5, are overplotted in red, green, 
and blue, respectively. 


The large velocity offset, high rotation rate, enhanced nitrogen 
(and possibly helium) all fit, at least qualitatively, with a scenario 
in which VFTS 358 is the former secondary of an interacting bi¬ 
nary that became unbound at the explosion of the primary star. 

6.3.2. VFTS 831 

VFTS 831 is remarkable in two regards. First as a luminous 13- 
type supergiant, likely to be the evolved descendent of a formerly 
more massive O-type star, and secondly because, from the defi¬ 
nition used in Table 4, it is located within SL639 (sal2"5 from 
the visual centre). With v T =210.4 ±6.2 km s -1 it also appears 
to be runaway compared to the other targets in SL639. A model- 
atmosphere analysis of the spectrum of VFTS 831 has been pre¬ 
sented by McEvoy et al. (2015), who classified it as an appar¬ 
ently single star. It is one of the cooler objects in the supergiant 
sample, with T e g = 14,000 K and log g = 2.1. 

Following the discussion in Section 6.2, VFTS 831 does 
not appear particularly remarkable in terms of stellar rota¬ 
tion. Indeed, it has a relatively low projected rotational veloc¬ 
ity and a substantial contribution from macroturbulent broad¬ 
ening (t) e sini, 0 rt) ~(40, 65)kms _1 , as found in most 13- 
type supergiants (see e.g. Markova et al. 2014; Simon-Diaz & 
Herrero 2014). From the intrinsically strong N II A3995 sin¬ 
glet transition, McEvoy et al. estimated a nitrogen abundance of 
12 ± log (N/H) = 7.94. Although large, this is consistent with the 
other single supergiants in the sample, which have abundances 
spanning approximately 6.8 to 8.1 dex. 

6.4. Single-lined binaries 

For completeness we also investigated the velocities of single- 
lined binaries for potential runaway systems. In the final col¬ 
umn of Table A.3 we include the observed range of velocities, 
6v r m ax, from Dunstall et al. (in prep.). Following the arguments 
laid out in Section A. 3 of the Appendix, we adopt an estimated 
error on the true systemic velocity compared to the single-epoch 


estimate of 0.15 x Sv rm ax- For stars with v T above (or below) the 
systemic velocity, we then subtract (or add) this estimated error, 
to give a corrected velocity estimate, v corr . Using the same 3er 
criteria discussed above we then looked for potential runaways. 

Only one single-lined system, VFTS 730, was found to be 
a potential runaway, with u corr =324.0 km s^ 1 (a 4.3 a outlier 
cf. the adopted systemic velocity). Its range of observed veloc¬ 
ities is relatively small (Sv rm a,x = 16.6 km s , thus only just 
above the detection threshold employed by Dunstall et al.) but 
it could be, for example, a long-period system with an eccentric 
orbit. Indeed, inspection of its velocity estimates as a function 
of epoch show them to be monotonically decreasing from 332 
to 316 km s -1 , suggesting it as a long-period system. Further in¬ 
vestigation of this (and other) binary systems will require com¬ 
prehensive spectroscopic monitoring to determine their centre- 
of-mass velocities and other orbital parameters. 

6.5. Origins of candidate runaways 

The production of runaways is thought to be from one of 
two channels, with evidence in support of both mechanisms 
(Hoogerwerf et al. 2000). In the dynamical-ejection scenario the 
runaway is ejected via gravitational interaction between single¬ 
binary or binary-binary systems in a dense star cluster (Poveda 
et al. 1967). In the binary-ejection scenario, the primary star ex¬ 
plodes as a supernova, leaving behind a neutron star or black 
hole. In most cases the compact object leaves the system de¬ 
pending on the magnitude and direction of the birth kick, and 
the companion is ejected with a velocity approximately equal to 
its orbital velocity (Blaauw 1961; Eldridge et al. 2011). 

The complex star-formation history of 30 Dor makes it hard 
to cleanly attribute our candidate runaways to either of the two 
scenarios. R136 is sufficiently young (1-2 Myr; de Koter et al. 
1998; Massey & Hunter 1998) that none of its members is 
thought to have undergone a supernova explosion. Thus, any 
runaways from R136 would have to be ejected via dynamical in¬ 
teractions. However, as noted above, the VFTS region contains 
several older subclusters and associations, including NGC 2060, 
which contains a supernova remnant (e.g. Chu et al. 1992), and 
Hodge 301 which was argued by Grebel & Chu (2000) as a likely 
site of tens of past supernovae from consideration of its mass 
function. 

One of the important diagnostics might prove to be the con¬ 
nection with the observed rotational velocities (Fig. 7), in which 
the apparently rapid rotators could have a binary origin. For ex¬ 
ample, the theoretical discussions by de Mink et al. (2013), and 
the suggestion that the high-i> e sin i tail of the distribution of 
the rotational velocities for the VFTS O-type stars could be ac¬ 
counted for by past binary interactions (Ramfrez-Agudelo et al. 
2013). Detailed dynamical analysis of the O-type stars in the sur¬ 
vey has been underway in parallel to the current work, and the 
broader question of the origins of runaway stars from 30 Dor is 
discussed in more depth by Sana et al. (in prep.). 

7. Summary 

We have presented comprehensive spectral classifications and 
radial-velocity estimates for the B-type stars observed in the 
VFTS. We now briefly summarise our main findings; 

- The members of the two older clusters, Hodge 301 and 

SF639, appear kinematically distinct from those in the 

younger clusters and general field population. 
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- The preliminary H-R diagram for stars in Hodge 301 and 
SL639 (Figure 5) suggests an age of ~10-15Myr for 
the population of SL639, slightly younger than that of 
Hodge 301 (15 ±5 Myr). 

- The systemic velocity for 273 single stars, excluding 
those in the two older clusters and candidate runaways, is 
*£ = 271.6 ± 12.2 (s.d.) km s _1 . 

- Employing a 3-cr velocity threshold, we identify nine sin¬ 
gle stars as candidate runaways (2.9% of the single stars 
with radial-velocity estimates). These are mostly unevolved, 
early B-type objects, but include a B5-type supergiant 
(VFTS 831). 

- There appears to a bimodal distribution of rotation rates in 
the candidate runaways. Five stars have v e sin i < 65 km s^ 1 , 
while the others have v e sin i > 345 km s _ 1 . 

- Excluding VFTS 831 (because of its evolved nature), there is 
a remarkably small probability (^0.001%) that the projected 
rotational velocities for our other eight candidate runaways 
could be randomly drawn from those for the full (unevolved) 
B-type sample. 

- VFTS 358 has the largest velocity offset from the sys¬ 
temic velocity of 30 Dor ((5y r =— 106.9 ± 16.2 km s _1 ) 
and appears to be rotating rapidly (v e sin i = 345 km s -1 ). 
Preliminary spectral analysis suggests evidence of signifi¬ 
cant nitrogen enrichment, perhaps indicative of past binary 
interaction. 

The classifications presented here will underpin the quan¬ 
titative analyses in progress on this substantial part of the 
VFTS data. For example, they will enable investigation of the 
mass/temperature dependence (if any) of the estimated v e sin i 
results from Dufton et al. (2013), particularly when also com¬ 
bined with the results for the O-type stars from Ramfrez- 
Agudelo et al. (2013). In the longer-term, the ongoing proper- 
motion studies using HST imaging (programme GO 12499, P.I. 
Fennon) will add another dimension to our studies of the com¬ 
plex dynamics in this region, by providing tangential velocities 
(see Sabbi et al. 2013). 

Our radial-velocity estimates should also lend themselves to 
detailed investigation of the structural properties of the VFTS 
sample (e.g. using the techniques discussed by Parker et al. 2014; 
Wright et al. 2014), particularly when combined with the results 
for the O-type stars from Sana et al. (2013); this should provide 
further insights into the formation and evolutionary history of 
the massive-star population in the region. 

Finally, there is an ongoing monitoring programme of 
FRAMES spectroscopy to characterise a subset of the O-type bi¬ 
naries (PI: Sana); a similar monitoring programme for the ~ 100 
binaries discussed here and by Dunstall et al. (in prep.) would 
be valuable to determine the full orbital parameters of the lower- 
mass population of binaries, while also providing robust esti¬ 
mates of centre-of-mass velocities to better explore their dynam¬ 
ics in the context of ejected runaways. 
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Appendix A: Ancillary material 

A.I. Spectral montages 

To illustrate the spectral sequences as a function of lumi¬ 
nosity class, examples of VFTS B-type spectra are shown in 
Figures A.I, A.2, and A.3. The principal lines used in classifica¬ 
tion, as summarised in Tables 1 and 2 are indicated in each fig¬ 
ure. Luminosity sequences at spectral types BO.2, B1 and B2.5 
are shown in Figures A.4, A.5 and A.6, respectively. For clarity, 
the spectra in all six figures have been smoothed and rebinned to 
an effective resolving power of R = 4 000. 

A.2. Detailed velocity estimates 

For completeness, we list the individual velocity estimates from 
each line, calculated using the rest wavelengths given in Table 3. 
Tables A. 1 and A.2 present the results for stars using the Set 1 
and Set 2 diagnostic lines, respectively; results for the binary 
stars (following the discussion outlined in Section 4.2) are given 
in Table A.3. 

A.2.1. Comparison with published velocity estimates 

The methods and line sets adopted to estimate radial veloci¬ 
ties (v r ) for the B-type stars differ to those employed by Sana 
et al. (2013) for the O-type objects from the survey. As a con¬ 
sistency check between the two studies, we analysed 45 of the 
apparently-single, late O-type stars (with 09.5 and 09.7 types) 
and compared our velocity estimates with those from Sana et al. 
We find a mean and standard deviation between the two methods 
of Av = 0.5 ± 6.2 km s -1 (where Ai> = v T — fsana)- 

The differential results (Av) for the 45 stars are shown in 
Figure A.7 (with the standard deviations added in quadrature). 
For the purposes of global analysis of the sample and identifica¬ 
tion of potential runaways, these results are in excellent agree¬ 
ment. As a further check of any systematic trends arising from 
the different methods, the lower panel of Figure A.7 shows the 
same differential velocities, but now as a function of estimated 
v e sin i (from Dufton et al. 2013); no obvious systematic differ¬ 
ence is present. 

A.3. Centre-of-mass velocities 

We have investigated the expected differences between our mean 
values and centre-of-mass velocities using a Monte Carlo tech¬ 
nique. We assumed a sinusoidal velocity curve (i.e. circular or¬ 
bit) and, adopting the observational sampling of Field A of the 
survey (Table A.I of Paper I), we calculated ranges of observed 
velocities (i.e. 5v r max ) and their resulting mean velocity. We as¬ 
sumed a random initial phase, and considered periods ranging 
from one day up to several years. 

These simulated mean velocities provide us with an estimate 
of the expected error compared to the genuine centre-of-mass ve¬ 
locity. For a given period, we can then average these differences 
over the range of assumed phases to find a typical error. The ratio 
of this typical error to <5u rmax is shown in Figure A.8. For bina¬ 
ries with periods of < 250 d, the expected error is approximately 
12-15 % of dv r Triax , with this increasing for longer periods, and 
then tending to a value of ~25 % for periods > 350 d (greater 
than the maximum time-sampling of the real data). 

These simplistic simulations (e.g. not taking into account po¬ 
tential eccentric orbits) have probably underestimated the un¬ 
certainties but, given the observational cadence of the survey. 
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Fig.A.7. Comparison of radial velocities (v r ) estimated using the 
methods in this paper with results from Sana et al. (2013), in which 
An = v r — vsana, and the plotted uncertainties are the standard devi¬ 
ations of both estimates added in quadrature. The upper panel shows 
the Av results in the sequence of VFTS identifiers (i.e. increasing right 
ascension), while the lower panel shows Av as a function of v e sin i. 
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Fig. A.8. Typical errors in the average radial velocities (calculated for 
the observational sampling of Field A of the VFTS ) for simulated bina¬ 
ries compared to their true centre-of-mass velocities (as a percentage of 
the expected range in radial velocities, Sv„nax, see text for details). 


they illustrate that the mean multi-epoch v T values should pro¬ 
vide reasonable estimates of the centre-of-mass velocities for 
short-period systems (<250d), which display relatively small 
peak-to-peak variations. 
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4000 4100 4200 4300 4400 4500 4600 4700 

Wavelength (A) 


Fig.A.l. Example B-type supergiants (Class la). The primary diagnostic lines for spectral classification are identified. In the spectrum of 
VFTS 525 these are: He II AA4542, 4686; Si III AA4553-68-74; Si IV AA4089. 4116; those identified in VFTS 269 are: He I AA4009, 4026, 
4121, 4144, 4388, 4471, 4713; Mg II A4481; Si II AA4128-32. The broad A4428 diffuse interstellar band is also evident in some sightlines. 
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Wavelength (A) 

Fig.A.4. Example luminosity sequence at B0.2. The identified lines are the same as those in Figure A.l; we note the greater intensity of the 
metallic lines with increasing luminosity (including the C ill + O II features at ~A4650). 



Wavelength (A) 

Fig.A.5. Example luminosity sequence at Bl. The identified He and Si lines are the same as those in Figure A.l. 
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C. J. Evans et al.: VFTS: Classifications & radial velocities for the B-type stars 



Wavelength (A) 

Fig.A.6. Example luminosity sequence at B2.5. The identified He and Si lines are the same as those in Figure A. 1 and the nebular lines have 
been truncated in VFTS 335 as indicated. 
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Table A.l. Line-by-line radial-velocity estimates for single stars using the Set 1 absorption lines. 
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Table A.2. Line-by-line radial-velocity estimates for single stars using the Set 2 absorption lines. 
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Table A.3. Line-by-line radial-velocity estimates for single-lined binaries using absorption lines in the LR02 observations. 


C. J. Evans et al.: VFTS: Classifications & radial velocities for the B-type stars 
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Table A.3. continued 


C. J. Evans et al.: VFTS: Classifications & radial velocities for the B-type stars 
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